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^ ■ Abstract 

■ Extreme times techniques, generally applied to non-equilibrium statistical mechanical processes, 

>>■ 



are also useful for a better understanding of financial markets. We present a detailed study on 
the mean first-passage time for the volatility of return time series. The empirical results extracted 
from daily data of major indices seem to follow the same law regardless the kind of index thus 



,1^ . suggesting an universal pattern. The empirical mean first-passage time to a certain level L is fairly 

\ different from that of the Wiener process showing a dissimilar behavior depending on whether L is 

O 

c/j ' higher or lower than the average volatility. All of this indicates a more complex dynamics in which 

Q ■ a reverting force drives volatility toward its mean value. We thus present the mean first-passage 

time expressions of the most common stochastic volatility models whose approach is comparable 
to the random diffusion description. We discuss asymptotic approximations and confront them to 
■ empirical results with a good agreement with the ExpOU model. 
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I. INTRODUCTION 



First-passage and extreme value problems are a crucial aspect of stochastic methods 
with a long tradition of applications to physics, biology, chemistry and engineering, all of 
them related to non-equilibrium processes. We thus can mention driven granular matter or 
polymers passing through pores J1J2I, Q], chemical reactions dynamics [^], reaction-diffusion 
systems js], coarsening systems jg, J fluctuatin g iri terfaces js, 9] or even nuclear fussion and 
light emission 10 1 among many others jlll . 12. Iisl. [l^ . In addition, many natural records 
need also a similar description as are, for instance, floods, very high temperatures and solar 
flares or earthquakes [isl ]. 

In studying extreme value statistics of a given time series one wants to learn about the 
distribution of the extreme events, that is, the maximum values of the signal within time 
intervals of fixed duration, and the statistical properties of their sequences. In hydrological 
engineering, for example, extreme value statistics are commonly applied to decide what 
building projects are required to protect riverside areas against typical floods that occur 
once in 100 years 15|. 

All this effort and knowledge have not substantially been introduced in the exploration 
of extreme events in financial markets and it has mostly remained inside physical sciences 
and engineering without any great spread outside them. Nevertheless, we believe that this 
perspective can be helpful by providing an alternative approach to extreme statistics that is 
different from that of the current mathematical finance IGj which can also result in a better 
control of the risk in financial markets. 

In the quantitative study of financial markets the volatility, originally defined as the stan- 
dard deviation of returns, plays an increasingly important role as a measure of risk. There 
are nowadays many financial products, such as options and other financial derivatives, which 
are specifically designed to cover investors from the risk associated with any market activity. 
These products are fundamentally based on the volatility, therefore, its knowledge turns out 
to be essential in any modern financial setting and, hence, in any financial modeling. 

One of the earliest financial models, the model of Einstein-Bachelier assumes that 
the volatility is constant being itself the diffusion coefficient of a Wiener process. However, 
this assumption is questioned by many empirical observations which are gathered together 
in the so-called "stylized facts" ISj. The overall conclusion is that the volatility is not 
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constant, it is not even a function of time but a random variable. Consequently the measure 
of volatility has become more difficult and questions like at what time the volatility reaches, 
for the first time, a determined value - which may or may not be extreme- are quite relevant. 

The main objective of this paper is therefore to study the mean first-passage time (MFPT) 
of the volatility process. We approach the problem both from analytical and empirical 
viewpoints. For one hand, we analyze the MFPT for daily data of major financial indices 
and observe that the MFPT curves of all indices follow an universal pattern when the 
volatility is scaled in a proper way. On the other hand, we obtain analytical expressions 
of the MFPT for a special class of two-dimensional diffusion models commonly known as 
stochastic volatility models in the quantitative finance literature jl9|. We next compare the 
analytical results with the empirical predictions which provides a test about the suitability 
of these analytical models. We incidentally note that these stochastic volatility models are 
analogous to the ones arising in the random diffusion framework [2^ and even to some 



multifractal models 



2l|. 



As mentioned, previous works on extreme times are, to our knowledge, scarce and 



mostly 



dealing with the return process but not with volatility. In our early works 22, l23|] we 
have analyzed the mean exit time of the return based on the continuous random walk 
technique and addressed exclusively to tick-by-tick data. Other examples studying the 



extreme time return statistics are given in Ref. [2^, |25| and specially in Refs. 26|, |27|] where 
the MET for the stock price is simulated using an stochastic volatility model as underlying 
process. And finally, there are also recent studies focused on the volatility data analyzing 
the interevent time statistics between spikes 28|, l29|] or the survival probability comparing 



the high frequency empirics with results from multifractal modeling |21l . |30| . 

We end this introductory section by pointing out that the analysis of extreme times is 
closely related to at least two challenging problems in mathematical finance which are of 



great practical interest: the American option pricing 31 



3211 and the issue of default times 



and credit risk 



33 



3J] . Both problems require the knowledge of hitting times, that is, first- 



passage times to certain thresholds. However, the typical mathematical approach there is 
quite different from the one we study here. 

The paper is organized as follows. In Sect. [Ill we outline the most usual stochastic 
volatility models. In Sect. IIIII we obtain the general expressions for the MFPT based on 
these models. In Sect. |IV]we analyze the averaged extreme time and examine its asymptotic 
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behavior. In Sect. |V] we estimate the empirical MFPT of several financial indices and 
compare it with the analytic expressions obtained in previous sections. Conclusions are 
drawn in Sect. |Vl]and some more technical details are left to Appendices. 



II. STOCHASTIC VOLATILITY MODELS 



The geometric Brownian motion (GBM) proposed by physicist Osborne in 1959 35| is, 
without any doubt, the most widely used model in finance. In this setting any speculative 
price S{t) is described through the following Langevin equation (in Ito sense) 

^ = ^dt + adW,{t), (1) 

where a is the volatility, assumed to be constant, fi is some deterministic drift indicating an 
eventual trend in the market, and Wi{t) is the Wiener process. However, and particularly 
after the 1987 crash, there is a compelling empirical evidence that the assumption of constant 



volatility is doubtful [18|, |36|, neither it is a deterministic function of time -as one might 
wonder on account of the non stationarity of financial data- but a random variable. The 
volatility is often related to the market activity [33']. In this way we are assuming that market 
activity is stochastic and governed by the random arrival of information to the markets. 

The hypothesis of a random volatility was originally suggested to explain the so-called 
"smile effect" appearing in the implied volatility of option prices 19|]. In the most general 
frame one therefore assumes that the volatility a is a given function of a random process 

Y{ty. 

ait) = a{Y{t)). (2) 

Most stochastic volatility (SV) models that have been proposed up till now suppose that 
Y{t) is also a diffusion process that may or may not be correlated with price, and different 
models mainly differ from each other in the way that a depends on Y{t). 

The usual starting point of the SV models is the GBM given by Eq. ([T]) with cr given by 
Eq. ([2]) and Y{t) being a diffusion process: 

dY{t) = F{Y)dt + G{Y)dW2{t). (3) 

In Eqs. ([T]) and ^ Wiit) {i = 1,2) are Wiener processes, that is, dWi{t) = ^i{t)dt, where 
S,i{t) are zero-mean Gaussian white noises with {^i{t)C,i{t')) = S{t — t') and cross correlation 



given by 

UtMt'))=rS{t-t') (4) 

(— 1 < r < 1). Incidentally we note that any SV model defined through Eqs. ([iMl) is, in 
fact, a two-dimensional diffusion process. 

The most common SV models in the literature are the following: 

(a) The Ornstein-Uhlenbeck (OU) model in which 



that is: 



See Refs. 



a = Y, F{Y) = -aiY -m), GiY) = k; 



dY{t) = -a{Y - m)dt + kdW2{t). 
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391] for further details. 



(b) The Cox-Ingersoll-Ross-Heston (CIR-Heston) model: 

a = VF, F{Y) = -a{Y -m^), G{Y) = k^] 



(5) 



then 



See Refs. 



dY{t) = -a{Y - m^)dt + kVYdW2it). 



4G 
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42[ for further details. 



(c) The exponential Ornstein-Uhlenbeck (ExpOU) model: 



a = me^, F{Y) = -aY, G{Y) = k; 



(6) 



hence 



See Refs. 



dY{t) = -aYdt + kdW2{t). 



43 



4J] for further details. 



(7) 



From the above equations we see that the volatility is also described by a one-dimensional 
diffusion process: 

da{t) = f{(r)dt + g{a)dW{t). (8) 
Thus, for the OU model a = Y and (see Eq. ([5])): 



da{t) = -a{a - m)dt + kdW{t). 



(9) 
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However, obtaining a differential equation for a{t) for CIR-Heston and ExpOU models is 
not direct, since in tfiese cases tlie volatility, a = cr{Y), is a nonlinear function of processes 



Y and the differentials of a and Y are connected through the Ito lemma 45|, |46|] : 



M>^) = ^^>^ + ^|^^^- (10) 



For the CIR-Heston model a = y/Y and from Eqs. ([6]) and ( ITOl) we get 

da{t) = -la(a- ^] +kdW{t), (11) 

where 

n = 77?,^ — 

4q; 



^2 

p = - — . (12) 



In the case of the ExpOU model a = me^ , and 

da{t) = -a \n{a/M)dt + kadW{t), (13) 

where 

M = me'='/2". (14) 

III. THE MEAN FIRST-PASSAGE TIME 

In this section and the next, we study the MFPT of the volatility process from an analyt- 
ical point of view. We postpone for a later section. Sect. |Vl the analysis of empirical mean 
first-passage times for several markets and their comparison with the analytical expressions 
obtained in Sects. IllltHVl 

Let us denote by Tx{a) the MFPT of the volatility process. That is, Tx^cr) represents the 
mean time one has to wait in order to observe the volatility, starting from a known value 
0", to reach for the first time a prescribed value A, which we often refer to as the "critical 
level". 

If we assume that the volatility is given by the diffusion process described in Eq. ([H]), 



then Tx{a) obeys the following differential equation [46 1 



2 da"^ da 
with an absorbing boundary condition at the critical level: 



ijV)^^ + /M?^ = -i (15) 



Tx(\) = 0. (16) 



Let us recall that the volatility should be a positive defined quantity. Hence, we have also 
to impose a refiection when it reaches the value a = 0. This is achieved by adding the 
refiecting boundary condition: 



Before proceeding ahead we note that the most general approach to the problem at hand 
would be obtaining the MFPT, T{S,a), of the two-dimensional process {S(t),a(t)) defined 
in Eqs. ([I])-©. After knowing T{S,a) -certainly a difficult mathematical task- we can get 
two different extreme times. Thus by averaging the volatility out of T{S,a) we have T{S), 
i.e., the MFPT for the price. On the other hand, averaging the price out of T{S,a) one 
gets the MFPT for the volatility T(cr), the latter being the main objective of the present 
work. Obviously T{a) is much easier to obtain from Eqs. (fT5l) - (fT71) than from this general 
and somewhat tortuous proceeding based on T{S,a). 

Let us return to the solution of the problem posed by Eqs. ffT^ - ffT7|) . This can be attained 
by elementary means with the result 



We shall now evaluate the expressions taken by the MFPT, Tx{cr), according to the SV 
model chosen. 

(a) The OU model. In this case (see Eq. (j9])) 




(18) 



where 




(19) 



fix) 



a{x — m), 



g{x) = k. 



Hence 



ip{x) 




2mx), 



(20) 



where 



a 



1/2 



u = 



k 



(21) 



Using Eq. with Eq. ([201), we get 




(22) 
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which in turn can be written as 



rz i'u(X—m) 

Txia) = — e^' [erf (z/m) + erf (x)] dx, (23) 

^ J u{<j—m) 



where 



2 2 
erf(x) = —= / dy 

VTi" Jo 

is the error function. 

(b) The CIR-Heston model. From Eq. (ITT]) we have 

/(x) = - ' 9{x) = k, 

and 

ijj[x) = -2z/^x^ + 4mz/^ \nx. (24) 
After substituting Eq. into Eq. (|T^ . some elementary manipulations yield 

Txia) = - x-^e'dx / y^-^e-^dy, (25) 



where 



P^lil + W)- (26) 



We use the following integral representation of the confluent hypergeometric function 47| 



F{a,c,x) = e-^x^-^ j\-yy^-'^-\l - xyf-^y, 

1 a 1 c — a 







and write Eq. ( l25i) in its final form 

Tx{^) = —5 F{l,l+P,x)dx. (27) 

(c) The ExpOU model. In this case (c/. Eq. ffT^ ) 

/(x) = — axln(x/M), (yf(x) = /ex. 

Consequently 

^{x) = -uHn^{x/M), (28) 

where z/ and M are given by Eqs. fl^ and ffT^ respectively. As before, substituting Eq. fl^ 
into Eq. (1181) and some elementary manipulations involving simple change of variables inside 
the integrals, result into 



2 rvW 

Tx{(j) = - e'^'dx e-y'dy, (29) 



a 



ri{a) 
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where 

ri{a) = ^ + u ln(a/m). (30) 
Note that using the error function defined above we can write Tx{a) up to a quadrature by 

Tx{(t) = — e"' [1 + erf (x)] dx. (31) 

« Jr){a) 

We finish this section reminding what is the MFPT for the simplest SV model. This 
is the case when the volatility is totally random without any reverting force driving the 
volatility to its normal level, that is 

da{t) = kdW{t), (32) 

where is a constant. We remark that, in contrast with the previous SV models, this 
dynamics has no temporal correlations at all, i.e., it has no memory. To compute its MFPT, 
we take again Eqs. ( fT8l) -( |T9l) . We insist in the fact the we are assuming a refiecting barrier for 
cr = in order to get a dynamics restricted between and A. This is in fact the reason why 
we obtain a finite MFPT since without a reflecting barrier which prevents a to be negative 



the MFPT would not exist 
a = 0, it is straightforward to get 

1 



46| . In such after including the reflecting barrier at 



This will be our benchmark solution in future sections. 



TM = i-A^'-^')- (33) 



IV. THE AVERAGED MFPT 

The expressions for T\{(t) developed in the previous section give us the mean time one has 
to wait until the volatility reaches a given level A startin g fr om its present value a. However, 



it is also of theoretical and practical interest [12, 122, l23|| the knowledge of the averaged 
MFPT in which the dependence on a has been averaged out. One might argue that this 
quantity has fewer applications to trading and investment but, as we will see later, for the 
current purposes of this paper this simpliflcation really helps to reach relevant conclusions 
based on real data. 

To obtain this average we have to choose a probability distribution for a. The simplest and 
most common assumption takes the present value of the volatility as uniformly distributed 



over the interval (0,A). For our purposes this choice is also very convenient because the 
average performed is independent of the SV model chosen, i.e., it is the same average for all 
models. We will thus test their abilities to reproduce real data under identical conditions. 
We thus define 

T{X) = jj^ Tx{a)da. (34) 

We can easily obtain T(A) for the SV models discussed above. This is done at once by 
substituting into Eq. (IH^ the expressions of T\{a) given by those models. For the OU, 
CIR-Heston and ExpOU models (c/. Eqs. f l2Hl) . fl?r|) and flHTl) ) this replacement, followed by 
an integration by parts, yields respectively 

T(A) = ^ / xe"'("-")' [erf (i/m) + eii{v{x - m))l dx, (35) 
aA Jo 

T(A) = -— / x^l^F{l, 1 + /3, x)dx, (36) 
apv\ Jq 



and 



T(A) = / +"/'^ [1 + erf (x)] dx, (37) 

aA 



— oo 



where, in writing the last equation we have used the definition of the parameter M given in 
Eq. fll4p . For the memoryless random volatility given by Eqs. (!H^ - fl33p we have 

T(A) = ^. (38) 

Equations (l35l) and (!36l) are the final expressions of the averaged MFPT for the OU 
and CIR-Heston models. The expression given by Eq. flHTl) . corresponding to the ExpOU 
model can be written in an alternative form which will show their usefulness both in the 
asymptotic and empirical analyses to be undertaken below. Thus, in Eq. ( 137|) we change 
the variable x —x and use the identity erfc(x) = 1 — erf(x) (erfc(a;) is the complementary 
error function). We have 

T(A) = r e"'-"/"erfc(a;)cix. 

o^A J-n{x) 

The complementary error function can be written in terms of the Kummer function U{a, c, x) 
as 1471 



1 _..2.. /I 1 2 



erfc(x) = — =e ^' ( 2' 2'^ 
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Therefore, 



which is our final expression of the averaged MFPT corresponding to the ExpOU model. 



A. Scaling the critical level 

Before proceeding with the asymptotic analysis of the above expressions of T(A), it is 
convenient to scale A so as to render it dimensionless. This will also be useful for the 
empirical analysis of the next section. 

In finance, it is rather relevant the knowledge of the so-called "volatility's normal level" , 
cTg, which is defined as the mean stationary value of the volatility process: 

/oo 
apst{(T)da, (40) 
-oo 

where Pst{c) is the stationary probability density function (pdf) of the volatility random 
process. The significance of as lies in the empirical fact that the volatility is, as time 
increases, mean reverting to its normal level 

We therefore scale the critical level A with the normal level, and define 

L^-. (41) 

From the analytical view, the dimensionless critical value L depends on the SV model we 
choose. The calculation of Us for different SV models is given in Appendix |X] with the result: 

f 

m OU model, 

= < 7/Z/V2 CIR-Heston model, (42) 
^gi/4i^^ ExpOU model. 



where 



and P is defined in Eq. 

Once we know the normal level, substituting A in terms of L into Eqs. fl35p . flHUj) and 
result in the following expressions for the averaged MFPT: 

(a) The OU model 
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T(L) = a^^^ f ^^m^'^Hx-if fgrf (m^y) + ed{mu{x - 1))1 dx. (44) 
aL Jq ^ 

(b) The CIR-Heston model 



1/2, 



= / 1 + (3, x)dx, (45) 



where u and 7 are respectively defined in Eqs. (I?T|1 and 
(c) The ExpOU model 



where 



C(L) = i- + ,.lnL. (47) 



Finally, we observe that the scaling of the critical level through the normal level ag given in 
Eq. ( HTl) cannot be performed on the memoryless model (j32ll . In this situation the volatility 
never reaches a stationary state and, hence, cTs is meaningless. If however we define a new 
critical level as L = A/A; (now L has dimension of time) we can write 

T(L) = h'. (48) 

B. Asymptotic analysis of the MFPT for small values of the critical level 

We will now obtain the asymptotic behavior of the averaged MFPT for small values of 
L. Note that the case L <^ 1 is equivalent to A ^ (Ts, that is, the critical level is a small 
fraction of the normal level. 

For the OU model T{L) is given by Eq. and one can easily see by means of a direct 
expansion that 

/ a;e'"''^'(""i)' [erf (mz/) + erf (mz/(x - l))]dx = -^L^ + O (L^) . (49) 
Jo SyTT 

Hence 

T(L) = ^L^ + 0(L3). (50) 



For the CIR-Heston Model, Eq. (I45l) . we use the expansion [47 1 



F(l,l+/?,x)=5^- 
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and get 



(51) 



The case of the ExpOU model requires a different approach than that of direct expansions. 
We first note that when L ^ the function C(-^) defined in Eq. (H71) tends to — oo. Hence, in 
this hmit the argument of the Kummer function f/(l/2, l/2,x^) appearing in the integrand 



of Eq. 



is exceedingly large. We can thus use the approximation 



47| 



X 



with the result 



where 



T{L) 



El [-C(^)/^] 



aL 



EAx) 



oo g-t 

— dt 
t 



is the exponential integral. Using the asymptotic approximation 



47| 



Ei(x) ^ 

X 



1 

1 + - 

X 



and taking into account (c/. Eq. (H7I) ^ 



— — + In L ~ In L, 



(L « 1) 



we finally obtain 



T{L) 



a In L 



+ 



In^L 



(52) 



Note that T{L) behaves, as L — > 0, in different ways depending on the SV model chosen. 
While in OU and CIR-Heston models (and in the memoryless model as well) the averaged 
MFPT grows quadratically with L, in the ExpOU model it grows logarithmically -an anal- 
ogous situation arises for large values of L, see below. We will return to this point in the 
following sections. 



C. Asymptotic analysis of the MFPT for large values of the critical level 

Let us now address the case L ^ 1 when the critical level is much greater than the normal 
level. In the Appendix [B] we show that an asymptotic representation of T(A) for large values 
of A and regardless the SV model chosen it is given by 

T(A) e-^^'\ (53) 

g^iXm)\/-2ip"{Xm) 
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Xr, 



where ip{x) is defined in Eq. (1191) . and Xm is the location of the maximum of ip{x). The 
values of Appendix [B]) 

m OU model, 

^ CIR-Heston model, (54) 
M ExpOU model, 

where p and M are given by Eqs. (|T2l) and (|T4l) respectively. 

For the SV models we are dealing with, Eq. fl53p up to the leading order of L yields (we 
use dimensionless units, cf. Eqs. fj^B - fHHj) ): 

(a) OU model 



where where z/ is defined in Eq. ( 12T|) . 
(b) CIR-Heston model 



T{L) ~ ^^^e^ " , (56) 



where 7 is defined in Eq. ( 14311 . 

(c) ExpOU model 



nL)-£yLe^''''\ (57) 



where we have assumed that — l/(4z/ ) + InL ~ InL for L ^ 1 



V. EMPIRICAL RESULTS 



We now present an empirical study of the MFPT for the daily volatility of major financial 
indices: (1) Dow- Jones Industrial Average (DJIA). (2) Standard and Poor's-500 (S&P-500). 
(3) German index DAX. (4) Japanese index NIKKEI. (5) American index NASDAQ. (6) 
British index FTSE-100. (7) Spanish index IBEX-35. (8) French index CAC-40 (see Table ID 
for more details). 

The empirical MFPT that we will obtain is the averaged MFPT and we will see its behav- 
ior in terms of the critical value. Let us recall that the uniform average has the advantatge 
over other more elaborated averages -as, for instance, averaging over the stationary distri- 
bution of the volatility- that the former is independent of the SV model chosen which allows 
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Financial Indices 


Period 


data points 


normal level (days ^/'^) 


DJIA 


1900-2004 


28,545 


7.1 X 10-3 


S&P-500 


1943-2003 


15,152 


6.2 X 10^3 


DAX 


1959-2003 


11,024 


8.4 X 10-3 


NIKKEI 


1970-2003 


8,359 


9.6 X 10-3 


NASDAQ 


1971-2004 


8,359 


7.8 X 10-3 


FTSE-100 


1984-2004 


5,191 


7.7 X 10-3 


IBEX-35 


1987-2004 


4,375 


9.6 X 10-3 


CAC-40 


1987-2003 


4,100 


10.2 X 10-3 



TABLE I: Empirical data used 



us to look at data on an identical footing when we confront the empirical observations with 
the predictions of any theoretical model. 

We work with the dimensionless level L defined in Eq. ( |4T1) so that we first need to know 
from data the volatility's normal level, cxs, of each market (see TableUfor the empirical values 
of (Ts). In this way we deal with critical levels that are proportional to the specific normal 
level of every market, thus unifying the magnitudes involved in the MFPT computation. 

We incidentally note that instead of the averaged MFPT T(A) we could have dealt with 
Tx(o"), that is, the MFPT starting from an specific value a and whose general expression is 
given in Eq. (ITSl) . However, as we have observed already in this case the statistics of data 
becomes less reliable. Moreover the casuistry in the data analysis become more complex 
and disentangling all their properties goes far away from our main objectives. 



A. Measuring the volatility 

From the time series of all the indices shown in Table [T] we have to extract first the daily 
volatility and then evaluate the averaged MFPT for many values of L. Before proceeding 
further we observe that, in fact, the volatility is never directly observed. In practice, one 
usually takes as an approximate measure of the instantaneous volatility (over the time step 
At assumed to be one day in our case) the quantity 

^ (58) 
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where AX{t) is the daily zero-mean return change defined as follows 



AX{t) 



ASjt) 



E 



AS{t) 



Sit) 



S{t) 



(59) 



where AS(t) = Sit + At) — S{t) are the daily price changes. The expected value appearing 
in this expression represents the conditional average of the relative price change knowing 
the current price S(t). 

We will now justify the measure of cr{t) given by Eq. (l58l) . Using Eq. ([T]) as the evolution 



equation of S(t), one can easily see that Eq. (13^ yields [48| 



AX{t) ~ a{t)AW{t). 



Consequently |AX(t)| ~ \a{t)\\AW{t)\; but \a{t)\ = a{t) and 



(60) 



\AW{t)\ = ^AWitf ~ At, 
where the last expression must be understood in mean-square sense since AW{tY — ^ At for 



sufficiently small At 



46| . Collecting results we obtain 



|AX(t)| ~ a{t)At 



which is Eq. ^ |49|. 

From the extreme-time problem point of view, we can check the soundness of this pro- 
cedure (i.e., identifying a{t) by |AX(t)|/At) comparing the results given by the exact an- 
alytical expressions of the MFPT for a obtained in Sec. [IV] with computer simulations of 
|AX(t)|. Unfortunately the evaluation of the MFPT T{L), for any critical lavel L, has to 
be performed by the numerical integration of Eqs. fj44l) - P6l) . The parameters we use for the 



numerical calculations are those given in the literature to reproduce the DJIA [39|, |42, |44| 
and they are summarized in Table HTl The methodology for evaluating the MFPT for |AX| 
from simulated time series is the same as that of the empirical data oulined few lines be- 
low. One should mention that the numerical computation of the integral is quite fast and 
straightforward except for the expOU model with a difficult numerical convergence from 
small values of L until the level raises the normal level of volatility. 

In Figure [T] we compare the two different results of the MFPT depending on the way 
we measure the volatility. The numerical form for a{t) may differ significantly from the 
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Parameters 


k 


a 


m 


r 


Ornstein-Uhlenbeck (OU) 


1.4 X 10^3 


5 X 10-2 


1.2 X 10-2 


-0.4 


Heston 


2.45 X 10-3 


4.5 X 10-2 


9.28 X 10-3 


-0.4 


Exponential OU (ExpOU) 


4.7 X 10-2 


1.82 X 10-3 


8 X 10-3 


-0.4 



TABLE II: The parameters, measured in daily units, for the SV models defined in Eqs. ([T])-©. 



OU 

Hesfon — 
ExpOU 
|dX| OU simulation * 
;dX .lation □ 
|dX| ExpOU simulation 


1 

^■4^ □ o 
^'■/ □ <> 






■ ^i"" X-::.-' 









0.001 0.01 0.1 1 10 



Proportion to normal level 

FIG. 1: Log-log representation of the MFPT for the OU, CIR-Heston and expOU models. Points 
are the MFPT results for |AX(t)| sample paths of 100,000 steps. The curves correspond to the 
numerical computation of the MFPT for a using the three integrals given by Eqs. (I44p - ()46p . 
Differences indeed appear but the qualitative behavior is preserved. 



simulation of |AX|. In all cases we observe that the simulated |AX| is noisier than a{t) as 
should have been expected because of the extra noise source /\W{t) appearing in Eq. ( l60l) . 
Thus the noise APU brings |AX| ~ cr|Aiy| to both larger and smaller values than those 
typically given by a alone. However, in the asymptotic limit the functional forms of |AX| 
and 0" are completely similar. 

Before going further, one should mention that the procedure just outlined to catch the 
true volatility a is not unique. There is, at least, another relatively simple way of extracting it 
from the time series of prices 4^] . This alternative method basically consists in dividing the 
empirical /S.X{t) -obtained through Eq. (l59l) - by a simulated Gaussian process replicating 
the Wiener increments APU(t). This, after using Eq. flUUj) . yields an empirical value for cr(t). 
We have shown in |4J| that such a "deconvolution procedure" works relatively well and it 
reproduces reliable values of the empirical volatility as long as the effects of memory and 
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cross-correlation are negligible (or, at least, that they do not affect the statistical analysis). 
However, in the analysis of extreme times the deconvolution method may destroy many 
subtleties of the MFPT curves. In this case, as we will see below, the market memory seems 
to really matter. 



B. The empirical MFPT 

Once we have a volatility time series constructed using Eq. (158|) . we can compute the mean 
first-passage time 5^. In Fig. [2] we present the log-log representation of the empirical T(L) 
in terms of L for the markets outlined in Table [B We observe that the MFPT behaves in a 
very similar way for all markets and universality is well-sustained. In the same plot we have 
added several curves. The solid one gives the MFPT for the volatility of the (memoryless) 
Brownian motion process: da = kdW{t). In this case, as we have shown in Eq. ( HHI) . 
the averaged MFPT is a quadratic function of L, T{L) = (2/3)L^, for all L > 0. Clearly, 
this behavior is not observed in real data. It seems to be therefore necessary to include a 
non uniform driving force in the drift that provides long correlations and clustering to the 
volatility dynamics. Moreover, we can fit two power laws of the form 

T{L) ~ L'^, (61) 

with a different exponent depending on whether L < 1 or L > 1 (c/Fig. [2]). 

As we have just remarked the quadratic law T{L) oc corresponds to a purely random 
volatility without memory. Moreover, for the standard GBM given by dX = adW, where 
X{t) is the zero-mean return and the volatility is constant, we see at once that the averaged 
MFPT for I AX I obeys the same law: 

T{L) oc (62) 

Therefore, following the discussion of the previous section about using |AX| as a measure 
of 0", the straight line with slope equals to 2 in Fig. [T] can be also understood as the case 
of the log-Brownian motion but for the returns increments itself. This could be seen as a 
benchmark and again stresses the need of including memory effects in any market model for 
getting the appropriate curves from levels greater and lower than the normal level. 

We also note that the range where statistics is reliable enough appears to be between 
0.01 and 10 times the normal level. The index with less statistics is the CAC-40 which is in 
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Proportion to the normal level 

FIG. 2: Log- log representation of the empirical MFPT, T{L), for the indices outlined in Table I, in 
terms of the dimensionless critical level L. Observe that the MFPT has a very similar behavior for 
all markets. The solid line represents the Wiener process for the volatility, that is, taking /(cr) = 
and g{a) = k, and for which T{L) = (2/3)L^ (as given by Eq. (|38|) ). Schematically we can also fit 
two very distinct regimes for L < 1 and L > 1 with T(L) = (2/3)L" whose exponent is respectively 
1.21 ±0.04 and 2.9 ±0.1. 
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FIG. 3: Log-log representation of the empirical MFPT, T{L), of the DJIA, along with the OU 



Heston and expOU MFPT's. Parameters of the models are the ones provided in the literature 



4|. 



fact the market with a lower data amount. The most reliable data statistics is that of the 
Dow Jones which allows us to look at data below 0.01 and far beyond 10 times its normal 
level. 

In the next figures, we compare the experimental result corresponding to the volatility of 



19 




0.1 1 10 
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FIG. 4: Log- log representation of the empirical MFPT, T(L), of the DJIA, along with several 
asymptotic adjustments corresponding to the theoretical SV models discussed in the text. 
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FIG. 5: Semi-log representation of the empirical MFPT, T'(L), of the DJIA, along with several 
asymptotic adjustments corresponding to the theoretical SV models discussed in the text. 



Asymptotics 


Fitting Region 


OU / CIR-Heston 


ExpOU 


L« 1 


< L < 0.4 
(25 points) 


a = 4.8 ± 0.2 day 


-6/lnL 
b = 0.59 ± 0.02 day 


L> 1 


1 < L < 10 
(91 points) 


cL eyip{uL'^) 
c = 1 day 
u = 0.9 


vLex.p{wln'^ L) 
f = 57 ± 1 day 
w = 0.17 ±0.02 



TABLE III: Fits for each SV model in the asymptotic regimes as given by the set of Eqs. (I50p - ()57p . 
We use a nonlinear least-squares algorithm except from the OU and CIR-Heston case for large L 
where it is not possible to provide a good fit and we just give some reasonable numbers. 
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Proportion to the normal level 



FIG. 6: Regular representation of the empirical MFPT of the DJIA for small values of the critical 
level L. The dotted line corresponds to the quadratic behavior T{L) ~ shown by OU and CIR- 
Heston models, Eqs. (j48p - ()5ip . Solid line corresponds to the logarithmic increase T{L) ~ 1/lnL 
of the ExpOU model, Eq. i^. 



Dow Jones 
expOU large L 
OU and Heston large L 




0123456789 10 
Proportion to the normal level 



FIG. 7: Regular representation of the empirical MFPT of the DJIA for a wide range of values 
L. The dotted line corresponds to the quadratic exponential growth, T(L) ~ e^^, shown by OU 
and CIR-Heston models, Eqs. ([5^ - dMI)- Solid line corresponds to the milder exponential growth 
T{L) ~ e^""^^ of the ExpOU model, Eq. §7^). 



the DJIA with the theoretical models and their asymptotic expressions. Figure [3] compares 



empirical data with simulations taking realistic parameters obtained in the literature 
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42 



44l |. We observe that the asymptotic regimes are better described for the expOU than the 



OU and Heston models. We can go a bit further but now taking the analytical expressions 
for a provided in previous sections. In Fig. HI we can see there that for both small and 
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large critical values, the theoretical model that seems to more accurately follow empirical 
data is the ExpOU model. The same conclusion is supported by the semi-log representation 
given in Fig. ([5]). For L ^ 1 this is enhanced in Fig. [6] where, in a regular plot, we see 
that the ExpOU model follows more closely and for longer values of L the empirical result. 
A similar situation is shown in Fig. [7] where the exponential growth of T(L) provided by 
OU and CIR-Heston models (Eqs. and (1561) ) deviates very quickly from the empirical 
MFPT, while the slower exponential growth of the ExpOU model (Eq. fl37|) ) seems to better 
fit the empirical result. In Table UTTl we show some details of the fitting procedure used to 
generate Figs. I1H71 

VI. SUMMARY AND CONCLUSIONS 

We have studied an aspect of the volatility process which is closely related to risk man- 
agement: its extreme times. The techniques employed have been mainly used in the context 
of physical sciences and engineering. In this way we have tried to broaden the field of 
applicability of the mean first-passage time (MFPT) to financial times series. 

We have estimated the empirical MFPT of several major financial indices. For all markets 
-from the American DJIA to the French CAC-40 and also for different periods of time (see 
Table [T])- the empirical MFPT follows an universal pattern as shown in Fig. [21 

These results sustain the assumed general random character of the volatility similarly 
to the random diffusion framework and they detect a very different behavior depending on 
whether the critical level is higher or lower than the average stationary volatility, i.e., the 
volatility's normal level. 

We have found that the MFPT versus the critical level L can be represented as a power 
law for each regime (see Eq. (16T1) and Fig. [2]). Compared to the Wiener process in which 
a = 2, when the critical level is below the volatility's normal level we get a ~ 1.2 which 
indicates a slower growth of the MFPT than that of the Wiener process. On the other hand, 
when the critical level is above the normal level the observed exponent is a ~ 2.9 and the 
growth is slower. All of this, has led us to look for theoretical models of the volatility which 
contain long correlations and are able to describe these distinct patterns depending on the 
average volatility. 

A second issue of our work has been devoted to the most common SV models whose 
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framework is analogous to that of the random diffusion approach. We have therefore chosen 
the OU, the CIR-Heston and the ExpOU models and obtained closed expressions (up to a 
quadrature) for the MFPT , Tx{(j), to a certain critical level A, being a the current value of 
the volatility (Eqs. ( l23l) . ( 127|) and ( l3Ti) ). By averaging out the current value of the volatility 
we have also obtained the averaged MFPT, T(L), where L is a dimensionless critical level 
representing the proportion to the normal level of volatility (Eqs. fl41l) - fl46l) ). 

Obviously different SV models furnish different expressions for the MFPT. However, the 
expressions obtained from OU and CIR-Heston models show a similar behavior while that 
of the ExpOU model is distinctive. This is clearly seen by asymptotic analysis. Thus, both 
OU and CIR-Heston models present, for small critical levels, a parabolic increase (Eqs. (H8ll - 

my- 

T{L) ~ L\ (L « 1); 
while for large values of L they show an "explosive" quadratic exponential growth (Eqs. fl3^ - 

T{L) ~ {L > 1). 

On the other hand the ExpOU model displays, for small critical levels, a logarithmic 
increase (Eq. (1521) ): 

T(L)~l/lnL, (L<1); 
and a milder exponential growth when L is large (Eq. (1571) ): 

T(L)~e'"'^, (L>1). 

Moreover, we have also fitted the above asymptotic approximations provided by SV 
models to the empirical data (see Figs. HHTj), with the overall conclusion that the ExpOU 
model better explains the experimental MFPT than OU and CIR-Heston models, specially 
for large values of the critical level. 
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APPENDIX A: THE NORMAL LEVEL OF VOLATILITY 



We know that the normal level of volatility, cXs, is defined as the mean stationary value 
of the volatility process (T(t): 



apst{cr)da, 



(Al) 



where Pst(o") is the stationary pdf of <j(t). For SV models the volatility process is determined 
by the one-dimensional diffusion given in Eq. ([H]), in this case the stationary distribution is 



explicitly given by 



46| 



Pst[cr) 



N 

2T^exp 



2 / Mda 



(A2) 



where iV is a normalization constant. For the SV models discussed in the main text, the 
stationary pdf reads 
(a) OU model: 



V 



(A3) 



where u is defined in Eq. f l21l) . 

(b) CIR-Heston model: 



Pst(cx) 



(A4) 



where f3 is given by Eq. 
(c) ExpOU model: 



(A5) 



where M is defined in Eq. ([T] 

Finally from Eqs. ( lAll) and Eqs. ( 1A3I) - ( lASp we immediately obtain the normal level of 
every SV model (c/. Eq. (pi): 



0"s = < 



m OU model, 

CIR-Heston model. 



(A6) 



me 



ExpOU model. 



APPENDIX B: ASYMPTOTIC APPROXIMATIONS 



In this appendix we will prove that a convenient asymptotic representation of the averaged 
MFPT, for large values of the critical level, is given by Eq. 
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The starting point of our derivation is the general expression of the MFPT given by 
Eq. ffT^ . We introduce this expression into the definition of the averaged MFPT given in 
Eq. (l34l) and then exchange the order in which the double integral is performed. We have 

= T / -::^dx I ye-^^y^dy, (Bl) 

which, after defining 



A Jo g'^ix) 



1 



h{x,X)^^ I'ye-^^y^y, (B2) 
9 [x) 



can be written as 



T(A) = ^ [ h{x,X)e^^''^dx. (B3) 
^ Jo 

Let us now see that, in all SV models studied in this paper, the function 



X{x) = e'^(^) 

reaches its maximum value at a point Xm which, for A sufficiently large, is inside the inte- 
gration interval (0, A) of Eq. fIBSp . Indeed, the extremes of x{x) coincide with extremes of 
ip{x), and from Eq. (|T9l) we see that the extreme points of iplx) are those in which the drift 
f{x) = vanishes. Therefore, for all SV models here analyzed, there is one and only one 
extreme given by 

m OU model, 

^/m CIR-Heston model, (B4) 
M ExpOU model. 

Note that for A large enough (in fact larger than the maximum value of m, \fm and M) 
Xm lies inside the interval (0, A). Moreover, Xm is a maximum since the second derivative 
■ip"{x) = 2f'{x)/g'^{x) < is negative for all models. 

The fact that the function x{x) = e^^^^ reaches its maximum value inside the integration 



interval (0, A) o: 
of the integral 



Eq. (1B3P allows us to apply Laplace's method for the approximate evaluation 



5l| . Expanding tp{x) around Xm 

ilj{x) = ^{Xm) + <(x„)xV2 + 0{X^ 



and substituting this into Eq. (1B3P we get 



2 '■^ 



T(A) ~ jh{xm, A)e'^("™) / e'^"("™)" /'cix. 
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FIG. 8: Approximate evaluation of /i(xm,A). 

but it is supposed that Q^'^'>P"(^m)/2 Qff quickly to zero. Hence, for sufficiently large values 
of A, we may write 

A foo rz 



Jo 2A/-V^"(xm)' 

whence 

T(A) ~ -^===h{x^, A)e^(^-). (B5) 

Now, in order to obtain an asymptotic approximation of T(A) as A — >■ oo, we have to 
discern the behavior of the quantity (c/Eq. ( ]B2p ): 

9 K^m) Jxm 

as A becomes large. One can easily show that the integrand of this equation, xe~'^^^\ is an 
increasing function of x for x > Xm (see Fig. [8]). We can, therefore, approximate the integral 
by the value of the area of the trapezoid depicted in Fig. [H] with the result 

h{xm. A) ~ 7^^^(A - Xm) [Ae-^W + x^e-^(^-)] , 

and for large values of A we can write 

Mx„,A)^— ^AV^W, (B6) 

which is the approximation sought for h{xm, A). Substituting this into Eq. (IBSp and taking 
into account that <^ |^(A)| (A 00) we finally prove Eq. 

r(A) e-^^'\ (B7) 

g'^{Xm)^/~2llj"{Xm) 
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which depends on the volatility and hence on the SV model chosen. We have nonetheless 
compared the results for the MFPT using these two ways of defining the return and we have 
not seen any significant difference. 

This procedure of extracting a from data furnishes only a "proxy" of the true value of the 
volatility. There might be more accurate analyses to obtain this value but any of these methods 
needs to assume an specific SV model which is one thing we want to avoid. In any case, the 
effort is far beyond the scope of this paper and we leave it for future future research. 
We first estimate the MFPT, Tx{cr), starting from a known value a and then raising the 
critical level A. For each critical level we take 200 boxes of the same width. For each box 
and each critical level we perform an average over the whole time series that allows us to get 
Tx{a). Afterwards for each critical level we average uniformly over all initial boxes and thus 
obtain the desired averaged MFPT given by Eq. (j34p . We perform all these computations for 
50 different levels both linear and logarithmically distributed between 0.001 until 100 times 
the normal level of volatility for the eight indices (see Table H]). The number of boxes and the 
number of levels have been chosen to be neither too big and nor too small. We get a reliable 
statistics having in all cases a small enough number of boxes to get a sufficient number of 
events inside each box and a large enough number of boxes for getting a good estimation of 
the averaged MFPT. 
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